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Abstract. We consider a modified version of the one-dimensional Hubbard model, the t\ —t2 Hubbard chain, 
which includes an additional next -nearest-neighbor hopping. It has been shown that at weak coupling this 
model has a Luttinger liquid phase or a spin liquid phase depending upon the ratio of t2 to ti . Additionally 
if the on-site interaction U is large enough, the ground state is fully polarized. Using exact diagonalization 
and the density-matrix renormalization group, we show that the transition to the ferromagnetic phase is 
either of first or second order depending on whether the Luttinger liquid or spin liquid is being destabilized. 
Since we work at T = 0, the second order transition is a quantum magnetic critical point. 

PACS. 71.10.Fd Lattice fermion models (Hubbard model, etc.) - 75.40.Mg Numerical simulation studies 
- 05.70.Jk Critical point phenomena 
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O . 1 Introduction there is no theory for this transition, and it should be in 



X 



the same universality class as the onset of ferromagnetism 

There has been interest in the theory of zero temperature in a Luttinger Uquid of itinerant electrons. ||] The criti- 

ferromagnetic transitions for a few years now. (D In such cal theory of this transition is the main remaining open 

transitions, it is the quantum fluctuations, rather than problem in the theory of phase transitions in quantum 

the thermal fluctuations, that govern the critical point, ferromagnetism. Q 

A theory for the onset of ferromagnetism in an unpolar- In this work, we will present a one-dimensional itin- 

ized itinerant system (Fermi liquid) for dimension d > 1 erant model which has a ferromagnetic quantum critical 

was proposed by Hertz who showed that the critical point. We study a modified version of the Hubbard model 

behavior should be mean field like. Precisely for d — I, by including a next-nearest-neighbor hopping in addition 
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Fig. 1. The ti — t2 Hubbard chain. 

to the nearest one. The model is no longer integrable, and 
to investigate it we used exact diagonalization and the 
powerful density matrix renormalization group (DMRG). 
1^ This model has previously been shown to have a para- 
magnetic to ferromagnetic transition as the on-site inter- 
action U is increased. Q Here, we show that the order of 
the transition is cither of first or second order depending 
on the parameters of the model. We will then focus on the 
second order transition. 

2 Model 

We consider the ti — t2 Hubbard chain (see Fig.|^) given 
by the Hamiltonian 

H = -hY^ (cl+iaCi<T + h.c^ 

-t2 {c\+2a'^i'y + ^-c-) + "»T"U- (1) 

The summation is over all L sites and spin a, and we will 
always take U positive. The sign of ti is arbitrary since a 
local gauge transformation, Cj — > e*'^^ Cj , maps the Hamil- 
tonian with ti negative onto the ti positive Hamiltonian. 
Therefore we set ti = \ without loss of generality, and 
measure all energies in units of ti. This Hamiltonian con- 
serves the number of particles N , the total spin S and its 
projection onto the quantization axis, 5*2. If a particle- hole 
transformation is applied to the system, the transforma- 
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Fig. 2. The non- interacting phase diagram. The dotted line 
indicates the boundary between the systems with two and four 
Fermi points. 

tion <2 is necessary to recover the original Hamil- 

tonian. We restrict ourselves to < iV < L and ^2 < 0, 
since it is in this region that a fully polarized ground state 
(with spin S — ^) has been found when U is large enough 
over a vast region of parameters. |6| We define Uc as the 
value of U above which the ground state is fully polarized. 

The weak-coupling phase diagram has two different re- 
gions which can be understood by looking at the U — Q 
phase diagram shown in Fig. ||. For not too large 1^2 1, the 
single-particle spectrum does not change much from that 
of the pure Hubbard. The Fermi surface has 2 points and 
we expect the model to be a Luttinger liquid Q away half- 
filling. When t2 < ^2^^! the Fermi surface has 4 points, 
and a weak-coupling treatment predicts that the system 
should be a spin liquid (no charge gap but a spin gap). 
[|j These two phases have been confirmed with DMRG 
calculations for [/ < Uc- @ 
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3 Order of the transition are needed in the Davidson procedure used in exact di- 



If the system is a Luttinger Hquid, it has gapless spin ex- 
citations with velocity Wg.. When it approaches the ferro- 
magnetic transition by increasing U , the velocity might 
smoothly go to zero leading to a second order transition. 
On the other hand, when the system has a spin gap, the 
transition must be of first order. To see this, following Ref. 
1^, a Hubbard-Stratanovitch transformation is performed 
on the interacting term of the Hamiltonian (|l|) . The intro- 
duced field m{q,Lu) can be seen as the order parameter, 



agonalization in order to obtain convergence (more than 
1000 H\t/j) multiplications). If the transition is first order, 
the ground state will jump from S* = to S* = S'max, and 
Eo{U) will have a kink at Uc, since the fully polarized state 
has no U dependence. On the other hand, if the transition 
is second order, the energy and spin will smoothly take 
on all values from to S'max as a function of U. In the 
thermodynamic limit, a second order transition requires 
that 



and the action of the system is developed in a power series u^Ua '^^ u^uj ou 

of m as i-e., the derivative of the ground-state energy is continuous 

through the transition, while it is discontinuous for a first 



S[m\ — + J dq duj [u2{q,CLi)m'^{q,uj) + UiUi^ + . . .] 

(2) 

where Sq is the non-interacting action, and u„ is the vertex 
of order n. The quadratic term of the effective action is 



order transition. In order to clarify this issue, we follow 
the lowest energy state with a particular spin S. However 
utilizing the quantum number in exact diagonalization 
is technically difficult, and so we follow a state of a par- 
^^2(9,^) = 1 — Ux{q,oj) (3) ticular S by diagonalizing the augmented Hamiltonian 



with x{<},^) the spin susceptibility. Since H' = H + AS^ (7) 

l^I^^I^Q-^^'^'^^ ~ ^ ^^-^ in different S^-subspaces with A > 0. For large enough A, 

for a spin-gapped system then 1°^°^* energy state within a given sector will have 

the minimum S value. ^ 

S[m] = So+ I dqdujm^{q,u;)[l + 0{q^,uj^)] . (5) 

J Results obtained with the Davidson algorithm Q for 

Since the coefficient of does not depend on U, the two different cases are shown in Fig. I In Fig. | a), for 

transition must be of first order. ^2 = —0.2, i = 12 and = 6, when the non-interacting 

Numerically we can study the order of the phase transi- Fermi surface has two points, we clearly see that the spin 

tion by calculating the ground-state energy Eo(U) around S of the ground state takes on all intermediate values as 

Uc with very high precision. Since there are many states U is increased. This is an indication that Eq. (^) will be 

with energy very close to i?o, a large number of iterations satisfied in the thermodynamic limit and that the transi- 
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Fig. 3. Ground-state energy as a function of U for different 
S'-subspaces. a) The system has L = 12, N — 6 and t2 — —0.2. 
b) The system has L = 14, iV = 8 and t2 = -0.8. Note that 
the sohd horizontal hue is the ground-state energy of the fully 
polarized state and does not depend on U. The state S = 3 is 
not shown in b) because it is higher in energy. 

tion is continuous. On the other hand, in Fig. |3|b) we see 
that the transition for a system with t2 = —0.8, L = 16 
and iV = 8, which is in the spin hquid phase at weak U, 
is from the S — Q state to the fuUy polarized one with 
= 4, indicating a first order transition. 
Another way of determining the order of the transition 
is to study the first derivative of the ground-state energy 
with respect to U directly. For the ti — 12 Hubbard chain. 
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Fig. 4. Double occupancy D as a function of renormalized U 
for two systems with L = 40 and = 20. The value of Uc is 
determined as in Ref. [5]. 



this is simply the double occupancy 
dEo 



dU 



(8) 



Here we perform DMRG calculations keeping up to 800 
states on lattices of up to 80 sites so that the maximum 
weight of the discarded density matrix eigenvalues is 10^^. 
Figure ^ shows the results as a function of renormalized 
U. We see that for t2 = —0.2 (full circles) the transition 
is continuous whereas when ^2 = —0.8 (open squares), it 
seems discontinuous. It is certainly the mixing of energet- 
ically close states near the transition which makes it seem 
smooth. 

4 Critical Exponents 

We will now focus on the second order transition. Here 
the relevant scaling field is not the temperature but the 
interaction 

g=\U-U,\. (9) 
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The order parameter m = (Sz) is coupled to the exter- leads to the identity 
nal field h. In addition to classical exponents one has to 
introduce a dynamical exponent z. The homogeneity hy- 
pothesis for the energy and the correlation function are 



z = d + 2-'i]. 



Thus for d — 1, one finds that v = (3 and z = 1 + Z. 



(19) 



(10) 



and 



rib-^/^g, b-yh, br, b'r) = 62-(<i+.+r,)/-(^^ ^) (H) 

with d the dimension, here d = \. This leads to the well- 
known definition of all critical exponents and to the fol- 
lowing identities 



a + 2/3 -f 7 = 2 
l = v{2~ 7j) 
a + (3{S + 1)^2 
{d + z)iy = 2 ~ a. 



(12) 
(13) 
(14) 
(15) 



with y related to S hy y ^ {d + z)S/ {1 + S). 

In addition, there is an identity derived by Sachdev. 
lllli Since the correlation function 



G(r,r) = (5, (r,T)5, (0,0)) 



(16) 



scales like 



G(6r,6V) =fe2-(rf+^+'?)G(r,r), 

the scaling dimension /i of Sz is 

d-2 + z + ri 
1^ = ■ 



(17) 



(18) 



However since Sz is a conserved charge density (it com- 
mutes with H), below the upper critical dimension its scal- 
ing dimension must be precisely d. Therefore, /i = d which 



The critical exponent 7 defined by 



(20) 



can be obtained in two different ways. The first is by tak- 
ing advantage of the fact that the system is a Luttinger 
liquid, for which the spin susceptibility is inversely pro- 
portional to the spin velocity. This leads for t2 = —0.2 
and n — 0.5 to a critical exponent of 7 = 2.0 ± 0.1. 
Another way is by adding to the Hamiltonian an external 
field h coupled to Sz 

H' = H + hSz, (21) 
where the susceptibility is then given by 

h^o ah 

The results for the same parameters as mentioned before 
are shown in Fig. ^. The critical exponent obtained by a 
least-square fit is 7 = 1.9 ±0.1 where the error comes from 
the fit. 

Another critical exponent can be obtained, namely a. 
Normally this exponent is defined by the divergence of the 
heat capacity. Since we deal with a quantum critical point 
this critical exponent is defined by 



zli?([/)=£;ferro-^o(C/)~5° 



(23) 



with a' = 2— a. Figure |^ shows the results for L — 40, t2 — 
—0.15 and n = 0.6. A mean square fit yields a' = 2.33 ± 
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Fig. 5. The spin susceptibility for L — 60, t2 — —0.2 and 
n — 0.5 as a function of U. The inset shows the points on a 
log-log scale fitted by the form Eq. (fed) with 7 = 1.9. 
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Fig. 6. The ground-state energy as a function of U for L = 40, 
t2 — —0.15 and n — 0.6. The dashed line is the energy of the 
fully polarized state. The inset shows the points on a log-log 
scale fitted by the form in Eq. (|2^) with a' — 2.33. 

0.05, where the error comes from the fit. A calculation for 
L = 80 yields the same exponent. 

Since these two exponents {a' and 7) only involve the 
evaluation of ground-state energies, they can be obtained 
with sufficient precision. From the identity between expo- 
nents, we get a small v which is consistent with previous 
calculation of correlation functions. B Table I shows the 



Table 1. Critical exponents for three different set of parame- 
ters. The errors are from the least-square fit. The critical ex- 
ponents f3 and z are calculated using the identities between 
exponents. 



t2 n 




a' 


7 


/3 


z 


-0.1 0.4 


15.1 


1.34 ±0.23 


1.15 ±0.16 


0.1 


12.5 


-0.2 0.5 


7.54 


1.87 ±0.03 


1.83 ±0.2 


0.08 


24 


-0.15 0.6 


16.3 


2.33 ±0.05 


2.17 ±0.06 


0.02 


100 



critical exponents obtained for three different sets of pa- 
rameters. We clearly see that they are not universal, and 
that /3 is small while z is large. This certainly is the sign 
of a crossover with the nearby first order transition (for 
which /? = 0), due to finite size effects. The fact that the 
first set of parameters, ^2 = 0.1 and n — 0.4, which are the 
farthest from the first order transition, gives the largest /3 
goes also in the direction of a crossover. 

5 Conclusion 

In conclusion, the t\ — t2 Hubbard chain, with negative 
<2 and filling less than one half, shows a transition from 
a spin liquid or a Luttinger liquid to a ferromagnet, de- 
pending on the ratio 12/ ti. This transition is of first order 
when we increase U from the spin liquid regime and sec- 
ond order when we increase U from the Luttinger liquid. 
This second order transition is characterized by a quan- 
tum critical point for which we can extract exponents. 
However the systems considered are relatively small and 
we see a crossover with the first order transition due to 
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finite size effects. More work has to be done by looking 
at even larger systems and also by trying to extract the 
dynamical exponent z directly. 
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